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Abstract 

Using soft-collinear effective theory we describe at leading order in l/m^ all the semi-inclusive 
hadronic B — > XM decays near the endpoint, where an energetic light meson M recoils against 
an inclusive jet X. We also include the decays involving n, rf mesons that receive additional 
contributions from gluonic operators. The predicted branching ratios and CP asymmetries depend 
on fewer hadronic parameters than the corresponding two-body B decays. This makes semi- 
inclusive hadronic B — > XM decays a powerful probe of the potential nonperturbative nature of 
charming penguins as well as a useful probe of new physics effects in electroweak flavor changing 
transitions. A comparison with B — » KX data from BaBar points to an enhanced charming 
penguin, albeit with large experimental errors. 
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I. INTRODUCTION 



Recently BaBar made the first measurement of semi-inclusive B — > KX branching ratios 
using fully reconstructed B decays [l| 

Br(B-/B° -> X-X) = (1961^+31) x 1(T 6 , 

Br(B-/B° - K°X) = (154+^+55) x 10 - 6j (1) 

where a lower cut on the kaon momentum p*(K) > 2.34 GeV in the B rest frame was 
imposed. This opens up the road for experimental explorations in hadronic semi-inclusive 
B decays, where for almost a decade the only observable probe has been Br (B — > tj'X), 
first determined by CLEO 0]. Averaging over the most recent measurements from BaBar 
[i| and CLEO 0] gives this branching ratio 

Br(£ -> 77%) = (420 ± 94) x 1(T 6 (2) 

for a lower cut on rf energy of E v r > 2.218 GeV. 

From the theoretical side semi-inclusive hadronic decays are very interesting since they 
are simpler, yet can still probe many of the questions that have been raised in the context 
of two-body B decays such as the perturbative and nonperturbative nature of charming 
penguins j5| and the search for new physics signals pj. Theoretical simplification occurs in 
the endpoint region, where the energy of the light meson M is relatively close to the maximal 
energy, so that the isolated energetic meson M and the inclusive collinear hadronic jet X go 
in opposite directions. Incidentally, this is also the part of phase space that is most readily 
probed experimentally. 

First predictions for the semi-inclusive hadronic decays B — > XM in the endpoint region 
were given in Refs. 0, @] using the soft-collinear effective theory (SCET) 0, 0, IE fiol (for 



earlier works on semi- inclusive decays using different theoretical approaches see [III, [121, ll3 , 



14], ll5|, [16|, [17], [18j). In the present work we go beyond Ref. [5j in several ways. First, because 
of the new experimental data on B — > KX branching ratios in Eq. (CQ), we are able to discuss 
the size of charming penguins and include it in perturbative predictions. Secondly, contrary 
to Ref. |i| in which only decays where the spectator quark is part of the inclusive jet were 
considered, we extend the discussion to all semi-inclusive decays, including decays to rj, rj ; . 
This is simplified by the fact that contributions where the spectator quark becomes part of 
the exclusive final state meson M are 1/?77| suppressed and can be neglected in our leading 
order calculations. These contributions are schematically shown in Fig. [lb to be compared 
with the leading-order contributions in Fig. [TK (additional gluonic contributions are present 
for decays into r/,r/'). 

This means that the nonperturbative parameters ( BM , (f M , connected to the B — > M 
form factors, do not enter in the leading order B — > XM predictions, making them simpler 



than the predictions for the corresponding hadronic two-body B decays 19j, |20j, |2l|, [22, 



23, 24, 25 



The presence of an inclusive collinear jet in the final state is described by a 
convolution of a nonperturbative shape function with a jet function. The latter arises in the 
matching of the full theory onto SCET! at the scale p 2 x ~ mbAq CD . At leading order this 
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FIG. 1: Time-ordered products of the effective weak operators for the decay widths: a) the leading- 
order contributions, b) a subset of subleading spectator interactions discussed in Section [TTJ The 
heavy quark fields are denoted by double lines, collinear quarks (gluons) by solid lines (overlaid 
with wiggly line) and soft quarks by dashed lines. The n hard-collinear quarks connecting the 
weak vertices and the n hard-collinear gluons [boosting the spectator in b)] carry p 2 ~ Ami, and 
are integrated out. The scaling of vertices is in A = y A/rrib, cf. Section ITU 

convolution is the same as in B — > decays, so that many hadronic uncertainties cancel 
by taking ratios. 

The paper is organized as follows. In Section [III we show the SCET power counting 
for the different possible decay contributions. This will allow us to include decays where 
the spectator could end up in the meson. We also discuss in this section the extra gluonic 
operators which contribute when the outgoing meson is an isosinglet meson. In Section HTT1 we 
briefly review the results of [5j and present the hard kernels for all semi-inclusive hadronic 
decays. In Section [TV] we discuss the production of r\ and rj mesons, where new gluonic 
operators are present at leading order in the power counting. In Section |V] we compare the 
predictions with data and then conclude in Section IVIl 



II. POWER COUNTING 

We work at leading order in l/my, as in Ref. At this order it is fairly easy to modify 
the results of Ref. j^] to include the semi-inclusive decays in which the spectator can go to 
either the jet or the light meson. In particular contributions where the spectator quark is 
boosted to become part of the exclusive final state meson M are 1 / m 2 suppressed and can 
be neglected. 

To show this we first explicitly power count different leading and subleading graphs in 
SCETi, where the expansion parameter is A = ^/A/m b . Order of the graph, X s , is given 



simply by power-counting the different vertices appearing in the graph 26 



5 = 4 + ^(A;-4)\4, (3) 

k 

where Vk is the number of vertices that scale as X k (the above equation already assumes that 
there are no purely ultrasoft diagrams). Now consider the contribution where the spectator 
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FIG. 2: Diagram that contributes to isosinglet meson production, which is of the same order as 
the leading diagram. 



ends up in the jet as shown in Fig. [TK For this graph Vq = 2 which gives a scaling 

5 jot = 4 + (6 - 4) x 2 = 8. (4) 

After SCETi — > SCETn matching there is an additional suppression of A for each external 
collinear line, which results in a final power counting of A 12 = (AqcdA 7 ^) 6 - This is the 
leading order term for semi-inclusive hadronic B decays. 

Next consider the contribution where the spectator ends up in the meson. A typical 
diagram is shown in Fig. [T]d with V4 = 2, V5 = 2, and Vg = 2, leading to 

Aspect = 4 + (4 - 4) x 2 + (5 - 4) x 2 + (6 - 4) x 2 = 10. (5) 

In addition to the usual A suppression for each external collinear line in the matching 
SCETi — > SCETn, however, this diagram is further suppressed due to the p± occurring 
in each of the A 4 collinear vertices. Lowering p± to the SCETn scaling gives an extra power 
of A for both vertices. This is exactly the same suppression that makes the soft-overlap and 
the hard-scattering contributions in heavy-to-light decays of the same order as discussed in 



Ref. [27j. The diagram in Fig. \T)p therefore scales as A 16 = (Aqcd/^&) 8 and thus is 1/m 2 . 
suppressed compared to the leading contribution and can be neglected in the leading order 
analysis. Other possible diagrams in which the spectator ends up in the final meson give 
the same suppression and can also be neglected. Note that this does not mean that all 
spectator interactions are 1/m 2 suppressed. In particular, annihilation contributions where 
the spectator quark annihilates with a collinear quark in the jet arise already at 1/rrib order 



as in B -> A7 [28 



The decays into isosinglet mesons 77, rf have additional contributions from gluonic oper- 
ators such as the one shown in Fig. [2} for which V5 = 4 and therefore 

5 iS o = 4+(5-4) x4 = 8, (6) 

which is the same power suppression as the diagram in Fig. After matching onto SCETn 
this diagram then contributes at leading power, (Aqcd/"^&) 6 - Thus, to analyze isosinglet 
meson production, we must include new contributions, complicating the analysis. Isosinglet 
meson production will be discussed in section [TVl 
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III. THE FORMALISM 



In this section we briefly review the results obtained in |5j while extending them to the 
full set of semi-inclusive decays. Additional contributions that arise for decays with 77 or 
7/ in the final state will be included in the next section. Barring those contributions, the 
decay rates of semi-inclusive B decays are obtained from the forward scattering amplitude 
of the time-ordered product of the heavy-to-light currents, as shown in Fig. [3j Because of 
disparate scales in the problem, a series of matchings on appropriate effective theories is 
performed. First, at the scale /1 ~ rrib the standard effective weak Hamiltonian in full QCD 
for hadronic B decays 29( is matched onto the effective Hamiltonian in SCETi by integrating 



out degrees of freedom of order m;, [19|, |20j, |2l|. In the next step SCETj is matched onto 



SCETn by integrating out the degrees of freedom with p 2 ~ m^A j^]. As a result, the jet 
function is obtained, the discontinuity of which contributes to the semi-inclusive hadronic 
B decay rates. This jet function is the same as in B — > X s ^ and will cancel once the 
ratio of decay rates is taken. The predictions for B — > XM branching ratios normalized 
to Br (B — > Xg'-f) and for direct CP asymmetries in B — » XM will thus depend only on 
perturbatively calculable hard kernels obtained from matching on SCETi at fi ~ rrif,, and 
on the remaining nonperturbative parameters - light cone distribution amplitudes (LCDA) 
and the parameters describing nonperturbative charming penguins. 

The hard kernels will depend on the Wilson coefficients Cf of the SCETi weak Hamiltonian 



that is at leading order (LO) in l/m^ given by JjJ, 120, \21 



H, 



p=u,c i=l 



(7) 



where ® denotes the convolution over collinear momenta fractions, while = V p bV* q is the 
CKM factor with q = s,d for AS" = 1,0 transitions. The Wilson coefficients Cf are shown at 
leading order in Appendix [A], and were calculated at NLO in 



a s [m bj 



first in Refs. 231, and 



then in Ref. 19|. In our notation, the NLO Wilson coefficients can be found in Appendix 




FIG. 3: Time-ordered products of Oi effective weak operators giving the decay widths through use 
of the optical theorem. The n collinear quarks connecting the weak vertices carry p 2 ~ Amj, and 
are integrated out. The gluonic contributions b) (with an additional mirror image not shown) and 
c) contribute only to B — * r]('X decays. 
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A of Ref . ^] . The sum is over four-quark operators 

1 = [uSPlXX] [lniPLUn] u , ^2,3 = [qJtP L Y*b v ] [v^Pl^]^ (8) 

q> q> 

and the gluonic operators are (the trace is over color indices) 



o 2g 



4tt 2 
9 2 m b 
An 2 



[qXynffn-Pig^P R Yib v ], 
[q n fiP L YX]n^B^] u ie ±llv 



(9) 



where the purely gluonic field is related to the (n, _L) component of the gluon field 
strength using the usual bracket prescription [20[ 

1 



n-P 



Wl[in-D n , iD£ ± )Wn 



(10) 



The operators 0\ g ^ g contribute only to the decays with rj, rf in the final state. We list the 
operator 2g for completeness, in order to have expressions valid to LO in l/m& but to all 
order in a s {m b ). When we discuss the phenomenology in section |Vj we work at 0(a s (nib))- 
At this order 2g has a vanishing matching coefficient 21] and thus does not contribute to 
the order that we are working. 

The summation over q' in Eq. (jSJ) includes u, d and s quarks and Pl,r = (1 =F 75)/2- 
The notation is the same as the one used in Thus [%^-Px%] u = \q n <5(w — |^)^P.r,<fo], 

while gauge-invariant n and n collinear quark fields q n = W^ffl and q n = already 
contain the collinear Wilson lines. The ultrasoft (usoft) Wilson line in the n direction, Y n , 
arises after the redefinition of the collinear fields to decouple collinear and usoft degrees of 



freedom 10 



This decoupling implies that the operators in Eqs. © and © factorize into currents 
Jc — (On ft^qQ and Jh = (q n f'PiXn^v) which do not exchange soft gluons. The matrix 
elements of n currents (M\Jc\0) are expressed in terms of light-cone distribution amplitudes 
(LCDA), while the time-ordered product of heavy-to-light currents, 



Mi 



m B 



d z e 



■ z (B\TJ}j(z)M0)\B), 



(11) 



leads to a convolution of shape, /(/+), and jet functions, Jp, |5| 



1 



7T 



ImT{E M ) = 2 



dl + f(l + ) 



-m b +2E M 



Im Jp ( nib — 2Em + 1+ + ie 



(12) 



= — S(Em, fJ-o), 



where l + = n • I is the soft momentum conjugate to the n ■ z = z_ spatial component. Using 
the optical theorem this is then related to the B — ► XM decay rate giving 



dY 



dE 



G 



-{B - XM) = ^mlx* M S(x M ,u. )\h^\ 2 + 



(13) 
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R~ — >• MY 


R° — » TtfY 


R° — > M Y~ 


rp(s) 

M,p 


A A UU 


7<"(*) _ 

/ \ v ' VY - 

ud 


A A us 


nV , nV 
L-l + '--4 


K ( * )0 X- 


K x dd 


A A ds 




a y — 




<P A S<S 


_i_ n _ _i_ ✓» 


VsX sa 




VsX sS 


nV \ n n 
1-4 + L-5 — L-6 


ujX su 


u>X°- 

sd 


uX% 


(CP + C 3 + 2C 5 + 2C 6 )/V2 


r lq X su 


^ X sd 


v q x% 


(C%-C 3 + 2C 5 -2C 6 )/V2 




- 0X %P 0X °sd 


n»X%,p°X% 


(Cf T C 3 )/V2 



TABLE I: Hard kernels 7gJ, with p = u,c for AS" = 1 semi- inclusive B~/B /B s -> IM de- 
cays. The NLO Wilson coefficients Cf are given in Appendix A of [&]. For the additional gluonic 
contributions to decays with r/ 9)S see section HVi 



5- - 


MA 




-> MX 




■* MX 


T (d) 



vr-A^^-Al n-X+p-X+ n~X+,p-X+ C* + C p 4 



*° X M>P° X M 


- 0x ^p° x °dd 


71 X ds> P X ds 


(C p 2 -CZtC 3 )/V2 


K^°X- U 


K (*)°x°, 

sd 


K^°X% 


ci 


" X du 


- x °dd 


" X ° d S 


(C% + CZ + C 3 + 2C 5 + 2C 6 )/V2 


^Q X du 


^ x dd 




{C p + Cl-C 3 + 2C 5 -2C 6 )/V2 


4 X dU 


^ x °dd 


4> x l 


C5 + Cq 


VsXJ a 


^ x °dd 


VsX% 


C5 — Cq 



TABLE II: Hard kernels T^' with p = u, c for AS = semi-inclusive B~/B /B s -> AM decays. 
The summation over p = u, c is implied. The NLO Wilson coefficients Cf are given in Appendix A 
of [H]. For the additional gluonic contributions to decays with 77^ see section HVl 

where Xm = ^Em/ttii, ~ 1. The ellipses represent nonperturbative charming penguin con- 
tributions given explicitly below, while h$ is the convolution of the hard kernel and the 
LCDA 




Here 4>m{u) is the light meson LCDA, f M the decay constant, X p 9> = V pb V* q the CKM 
elements, while the perturbatively calculable hard kernels are given in Tables [I] and [III 
for AS — 1, (q — s, d), respectively. 

The nonperturbative function S denoting the convolution of shape and jet functions is 
exactly the same as the one appearing in the prediction for the B — > X s 7 rate in the endpoint 
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a)- 





M 
cc 



FIG. 4: Charming penguin contributions: the contributions from a) and b) are proportional to p\ 
and pM* respectively, while the contribution from c) gives . The blobs represent nonperturbative 
charming penguins 



region at LO in l/m&. In the ratio with T{B — > X^) it thus cancels out, giving 
dT{B - MI)/^ _ 2.3 (|4T + 2^e[A^c cc p-(4?)1 + ^c^V") 



dT{B -> X s ~/)/dE~ 



ami 



A t W C 7 (cf + l/2c' 



(15) 



12/ 



where one sets EL = The SCET Wilson coefficients are cS ff = 1, 



c 



eff 
12 



at LO with 



'7 — -^M- 

the NLO calculated in [8(, while C 7 is given e.g. in Eq. (13) of 281 . 

The coefficients p^ c and parametrize possible nonperturbative charming penguin 
contributions. 1 They depend on the valence quark structure of both M and the jet X, 
but we suppress this dependence in the notation. They are zero if the charming penguin 
contributions are purely perturbative. However, the uncalculated higher-order perturbative 
pieces can mimick their effect, making them differ slightly from zero. The complex parameter 
p^. describes the interference of the nonperturbative charming penguin with the perturbative 
hard kernels, shown in Fig. H^,. The positive real parameter in Eq. (fl5|) on the other 
hand describes the square of the nonperturbative charming penguin contributions shown 
in Fig. Hb. If hard kernels dominate the amplitudes, the term with p^ in Eq. (TI5T) is 
subleading, while the term is even more suppressed and can be neglected as was done 
in jEf . It should be kept, however, if nonperturbative charming penguins are sizable. Since 
the present data are inconclusive we keep both terms in Eq. ( ITS"]) . When estimating the size 
of the non-perturbative contributions, as a rule of thumb we will take {p^.) 2 ~ Pir . Further 
information on the structure of p^,V^ can be obtained in the m c — * oo limit [5J. Finally, 
the coefficient c cc multiplying the nonperturbative charming penguin parameters in Eq. ( f!5l) 
is equal to the coefficient of C\ in Tables fl] and [IT] (i.e., it is c cc = 1 for B~ — > 7i~X® u and 
c cc = — l/v2 for B~ — > tt°Xj u so that p^f in both cases equals p^,). 

In the phenomenological analysis of our results in Section |V] we give numerical estimates 
for direct CP asymmetries 



Acp(B -> IM) 



dr(5 -> XM)/dE M -dT{B -> IM) /dE 



dT(fi -> XM)/dE M + dT(B -> XM)/dE M 



(16) 



In [6[ these were p, 



AT 



/m/cc and 7> c ^ = /f f ^ cc . 
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and CP averaged branching ratios. In terms of hard kernels and nonperturbative charming 
penguin parameters the direct CP asymmetry is 

_ -21m(X^X^*) r 

A C p(B -> XM) = _ K J lm\T(P + c cc p^Y] , (17) 

V 1 \A(B^XM)\ 2 1 V ycc) J ' V ; 

and the CP averaged decay width normalized to B — > X s 7 is 

dV C p(B -»• MX)/dE M 2tt 3 |A(5^XM)| 2 



(18) 



(19) 



rfr(5 -> X sl )/dE, ami {X^C^cf + l/2cff)| 2 ' 

where 

-> XM)| 2 =|Al 9) | 2 |T| 2 + 2ne(X^X^*)TZe[T(P + c cc p™)*] 
+ |AW| 2 0P| 2 + 2^e(c cc ^P*)+c c ^]. 
Above a shorthand notation for the perturbative "tree" and "penguin" contributions 

T = [ duf M <f> M {u)T$ u (u), P= f duf M cf> M (u)T$ c (u), (20) 
Jo Jo 

has been used, dropping in the notation the dependence on B — > XM. 

IV. THE DECAYS INVOLVING r/, rf 

In order to describe B — > r/'^X decays several modifications of the results in the previous 
section are needed: (i) 77 — 77' mixing needs to be taken into account and (ii) there are 
additional contributions from gluonic operators 0\ g ^ g as shown in Figs. EJ [5] [operators 
@ig,2g are defined in Eq. (J§J)]. To describe matrix elements involving gluonic operators we 



introduce the gluonic LCDA 30j, |31 



ie ±flu {P(p)\Tr[B^B^] u \0) = Z -^f P ^ P (u), (21) 

where the isosinglet decay constant is the same one that appears in the matrix elements of 
quark bilinears 

(P(p)m^T lMn ] u \0) = ^iEf/^/iu), (22) 

with T 8 = X$/V2, T\ = l/y/3 diagonal 3x3 matrices in u,d,s flavor space. The flavor 
singlet LCDA <p P (u) mixes with the gluonic LCDA $p(u) under RG running 32J, while 



!>p(w) does not. For future reference we also quote explicitly = ^2/3f Vq , /* s = f Vs /V^, 
while fp — for other pseudoscalars that do not have flavor singlet component. Here f Vq 
and f Va are the decay constants corresponding to qq = (uu + dd) / \/2 and ss axial currents 
respectively (and are equal in the SU(3) limit, cf. also (1251) below). 

The parametrizations of matrix elements (1211) . (I22p do not involve any assumptions; 
however they are too general for the limited amount of data available at present. To reduce 
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the number of unknowns we use the FKS mixing scheme 33j to describe 77 — 77' mixing in 
which the mass eigenstates 77, 77' are related to the flavor basis through 

T] = Tjg COS If — Tj s SU1 <f, f] = Tj q Sm if + Tj s COS <f, (23) 

with <p = (39.3 ± 1.0)° and rj q = (rj u + r]d)/\/2. The working assumptions of the FKS scheme 
are that LCDA do not depend on the meson so that 

<pf{u) = ^(u), & P {u) = (24) 

and that OZI suppression is effective. This last requirement is most transparent in the 
r) q , T) s , g basis instead of the 1, 8, g basis used above. In it we have 

(Vq(p)\^([Un^l5Un]u + [4 ^75^n]«) 1 0) = ~2lE f^^u) , 

(25) 

(Vs(p)\[Sn^l5Sn}u\0) = ~2iEf Va (j> Vs (u) , 

with <f) Vs (u) = [2<p 8 (u) + <p l {u)}/3 and <f) Vq {u) = [2(j>\u) + <p 8 {u)]/3. The OZI suppressed 
matrix elements on the other hand are 

(Vs{p)\ J m{[Ufi$ r )bUn]u+ [dn$lbdn]u)\$) = ~2iE f Vs <f) opp (u) , 

v (26) 

(Vq(p)\[Snftl5S n ]u\0) = ~2%E f nq (f) opp (u) , 

where 4> opp (u) = V2[<p\u)-<f) s (u)]/3 and is negligible as long as 1 (m) — <p 8 (u). This relation 
is exact for asymptotic forms of LCDA, while it can only be approximate for physical values 
of \i since 1 (m) and <fr 8 (u) have different RG runnings, spoiling the relation for smaller 
values of /x. Phenomenologically, however, for \i above 1 GeV the relation is well obeyed at 
a percent level 3(3 . 

We are now ready to write down the results for contributions to B — > r](j]')X decays 
corresponding to Figs. [3] and HI These are described by Eq. (fl5l) but with and charming 
penguin parameters as given below. Utilizing the FKS scheme with Eq. (125p and setting the 
OZI suppressed matrix elements Eq. (126]) to zero, the functions in Eq. (Tl5|) are 

/#> = cos tpf^ ® X^T^ p - sin ® + h% (27) 

hf = sin^ f Vq cj> Vq ® AWrW, + ccstpf^ ® AWrW, + h 9 v „ (28) 

where ® denotes a convolution, while T^ 3:P are listed in Tables H]HT1 and the sum over p = u,c 
is understood. The gluonic contributions h 9 M coming from the Oi g operator insertions as 
shown in Figs. [3b, [3b, are zero to NLO in a s (m b ), i.e. to the order we are working. Explicitly, 
they are 

K„ = ~ x t ] ^fLK>) ® c ^ ( 29 ) 

for AS = 1 decays B q i — > T) s>q X sq >, while for AS = decays B q i — > r] s ^X dq -i we need to 
replace Xf —> Xf . The expressions for 77,77' final states are easily obtained using ( J23l) . As 
already stated, C 29 (m) = at NLO in a s {m^). 
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FIG. 5: The contributions of 0\ g operator to B — » MX q decays. The mirror image of Diagram a) 

i) 



is not shown as well as not the diagrams with C^\t replaced by C^\g-> °f- Fig. El 



For the charming penguin parameters in Eq. (|15[) . we make the following replacements 
for AS" = 1 transitions 

CccPlc -»• - sin <p p v c 3 c , c 2 cc V v cc -> sin 2 yV£ , 

(30) 

CccPcc -> COS V? p£ , C^P^ -> COS 2 <£>7^ , 



and for AS* = transitions 



^ccfcc it: Vcci ^cc' cc r> ' cc 



v/2 



CC ) ^CC ' CC o cc J 

(31) 



We next move to the contributions from the 0\ g operator, shown in Fig. [5l These 
contributions lead to a modified factorization between n and n degrees of freedom because 
of the additional soft gluon that is emitted at a light-like separation from the weak vertex in 
the n direction. Relegating the details to Appendix [Bj we quote here only the main results 
starting with the T products in Fig. [6] 

g H M = {MX\ijd A xT{O l9 {Q),C^ {cg) {x)}\B). (32) 

Then + Q cg describes the contribution of the operator O g to the decay into a color-singlet 
state. 

In the SCETj — > SCETn matching the intermediate hard-collinear gluon carrying p 2 ~ 
Amb is integrated out leading respectively to the jet functions Ji(u, and J g (u,k-) for 
the diagrams in Fig. [B^ and [Hb- Following the usual redefinition of fields the n collinear 



quark (gluon) lines decouple from the soft fields [lOj and lead to quark (gluon) LCDA after 
taking the matrix element. The soft Wilson lines Y n arising from the field redefinition and 
the soft gluon field emitted at x are taken to be part of the heavy-to-light current 

J H (0,x + ) = q n Y^Y n (0)^gM(x + )P R Y^b v (0), (33) 
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o lg 



FIG. 6: Time-ordered products between 0\ g and the subleading interaction terms in SCET La- 
grangian. The soft gluon (curly line) is absorbed in a definition of new B meson shape function 
f g . The intermediate hard-collinear gluon has offshellness ~ m^A and is integrated out. At LO in 
a s (\/ Anib) Diagram b) does not contribute because of the antisymmetry of gluonic LCDA. 



where we define 



(34) 



Note that the heavy-to-light current Jh depends on x + , since the soft gluon gets emitted 
away from the weak vertex. The 0\ g contribution to the B — > MX matrix element then 
takes a form of a convolution over both the soft momentum k- (conjugate to the position 
x + ) and the hard momentum fraction upM, 



du 



dk— dx^ 

47T 



■e- ik - x +/ 2 F M (k„,u)(X\j H (0,x + )\B). 



(35) 



The hard-collinear kernel multiplied by the LCDA is explicitly 



F (k-, u) 



47T 



/m0m(w) 



Ji(u,k-) Ji(u,—kJ)* 



u 



u 



+ y/C;&& M (u)J g (u,k. 



(36) 

and is in general both a function of the n momentum fraction u and the soft momentum 
At tree level, however, it is a simple product of functions that depend only on u and only 
on k_, 

F M (k„,u) 



I tree 



asmb t a i \( l , l \i n 
— ; — Jm<Pm{u){- + -)Ji{k- 

47T U U 



(37) 



tree 



since Ji(u, k-.)\t Tee = Ji{k-)\tree — l/(Nk-) and J g (u, k^)\ tvec = l/k- are independent of 
u. Furthermore, since $p(w) is antisymmetric, $p(tt) = — $p(u), the contribution from 
J g (u, k_) vanishes at this order. 

Therefore at least at leading order in a^VAmj) the additional gluonic contributions can 
be cast in the same form as the expressions for the decay widths for nonisosinglet final states, 
Eq. (TI3"]) . by moving the ^-dependent part into the definition of hard kernels, while including 
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the dependence on in the definition of the modified heavy-to-light current 
.7(0) = f ^±e- ik -^ 2 J H (0,x + ) J^k. 



I tree 



(38) 



Because of this simplification we will show in this section only the result for B — > XM 
decay width at leading order in a s (-\/Am b ), while the result valid to all orders in a s (\/Tvrnb) 
is given in Appendix [B] 

As in Section II I II we relate the B — > XM decay width to the time-ordered product of 
heavy currents using the optical theorem. We denote the T-product coming from a single 
0\ g insertion, shown in Fig. [5^, as 



M 



I 



<fz{B\TJ]j{z)J{Q)\B), 



(39) 



where Jh{z) = e 1 ^ mbV " > ' z {q n ^PLY^b v ){z). For the contribution coming from two insertions 
of Gig, shown in Fig. \5jp, we similarly define 



f gg (E M ) = — [ d 4 z{B\Tj\z)J(S))\B). 
m b J 



(40) 



In evaluating the time ordered product T g (E M ) we use the fact that n collinear quark 
fields do not exchange any soft gluons at LO in l/m b with the other fields in the T-product. 
Using the standard definition of the n-collinear jet function, 



(Q\Tq n {z)Um) = ^ z +) 52 ^ I ^ ± e- iK + z -/ 2 J P (K + + ie), 
2 J 2tt 



(41) 



and a shape function that, unlike f(l+) in Eq. (fl2l) . depends on two soft momenta because 
of the additional nonlocal structure present in J due to ^4 US , 



J dl + e il+z ~/ 2 J dr.e ir - x+ ' 2 f g {l + ,r-) = 

(B v \b v Y4z^)Y r lY n (0)^gAi(x + )P R Yib v (0)\B v ), 

the discontinuity of T 9 (Em) with respect to the intermediate states is given by 

1 



(42) 



Disc. T 9 (Em) = 2 j dl+drJIm 
2 

= — S g {E M , Ho). 



71 



-Jp(l + +m b — 2E M + ie) 



/ ff (Z + ,r_)Ji(r_)| tre e 



(43) 



Note that the new shape function f g (l+, r_) can be in general complex, hence S g can also be 
complex. However for decays rates, the complex conjugate should be added and the decay 
rates becomes real. 

Using the optical theorem the contribution to the decay width from one insertion of Gi g 
at LO in a s {y/Km b ), with the contribution of the mirror image of Fig. [5k., is therefore 

dT 



dE 



M' 9 



(~i2 

2 2 fl 

— ~ A m b x MJM 
tree 47T 



(44) 



x 2UeUct>M®^ ) T^ p ) 



M ® 4 q) Cl„^- ( - + - 



An \u 



1 



it 



Sg(E M , fl>0) 
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where C\ g is the Wilson coefficient of 0\ gy which is 1 at leading order. This expression is 
similar to Eq. (1131) . The hard kernels (convoluted with LCDA) in the curly brackets do not 
depend on soft momenta k_ and similarly the nonperturbative "shape" function S g does not 
depend on the large momenta fractions u. This factorization of the u and k_ dependence 
is a consequence of a special form of F M at leading order in a s (\/ Am b ), Eq. ( 1371) . and may 
not be present at higher orders, cf. Eq. ( |36l) . 

The contribution T gg , Eq. (HUj) . coming from two insertions of 0\ g leads to a shape 
function that depends on three soft momenta because of two soft gluon insertions 



(45) 



dl + e ll+z - /2 J dr_e ir ~ x+/2 J ds_e ls - y+/2 f gg {l+,r_, s_) = 
(5„|&^_)^(y + )^V^ 
This gives a new nonperturbative "shape" function by taking the discontinuity 



Disc. T g g(E M ) = 2 / dl+dr^dsJIm 



Jp(l+ + m b - 2E M + ie) 

7T 



(46) 



x J 1 (r_)J 1 (s_)*f gg (l + ,r_,s_) = — Sgg(E M ,fM}), 

mb 



which then enters the prediction for the contribution of double 0\ g insertion to the decay 
width 



dY 



AE 



M / 99 



G 



^mlx M fMS gg {E M ,fio) K q) Ci g I ducf)M{u)—-[ - + _ 

tree Z7T / 47T \U U 



« s / 1 1 
) — I h - 



(47) 



Note that f gg and S gg are real in contrast to f g and S g . 

To recapitulate, the prediction for the B — > T](j]')X decay widths at LO in l/m^ is given 

by 

dT / dY \ / dT \ / dT \ 

VdEA// Eq.GS) \dE M )g \dE M hg 
with the first term interpreted according to the replacement rules given explicitly in 



Eqs. (B7I)- (|3T1) . while the last two terms are given at leading order in a s (-y/Am&) in Eqs. (T45 
and f|47|) and to all orders in Appendix [B] 



V. PHENOMENOLOGY 

We are now ready to use the expressions for CP averaged branching ratios and direct 
CP asymmetries derived in the previous two sections for quantitative analysis. We split the 
discussion into two parts, first focusing on the decays to nonisosinglet pseudoscalar and to 
vector final states and then moving to the predictions for the B — > r)X, rj'X decays. 

While the first measurements of B — > MX decays have become available, one still lacks 
enough experimental information to determine nonperturbative charming penguin parame- 
ters from data (or to show decisively that they are small and compatible with zero). There- 
fore we collect in Tables II III and IIVI only purely perturbative predictions for B — > XM 
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MX 


Br(B —> MX)/Bi(B —> 


X sl) 


A C p 


K ~ X uu 


0.17 ±0.09 ±0.06 




0.30 ±0.16 ±0.01 




0.20 ± 0.11 ± 0.06 




(9.7 ±4.8 ±0.6) x 10~ 3 




(1.0 ±0.6 ±0.2) x 10" 


o 




K * ' X uu 


0.28 ±0.16 ±0.06 




0.32 ±0.16 ±0.02 




0.34 ±0.19 ±0.07 




(8.4 ± 4.6 ± 1.9) x 10~ 3 




0.22 ±0.13 ±0.03 




(8.9 ± 5.0 ± 1.6) x 10~ 3 


UjX SU 


(2.8 ±3.3 ±0.7) x 10' 


-3 


0.49 ± 0.24 ± 0.33 


P° X su 


(2.4 ± 1.4 ±0.5) x 10" 


-2 





TABLE III: Predictions for decay rates and direct CP asymmetries for charged B~ — » MX AS = 1 
semi- inclusive hadronic decays, which are the same as for corresponding B — > MX, B s — > MX 
given in Table [H The first errors are an estimate of the 1/mb corrections, while the second errors 
are due to errors on the Gegenbauer coefficients in the expansion of the LCDA. 



decay rates and direct CP asymmetries using Eqs. ( 1151) and ( jl6l) . setting the nonpertur- 
bative charming penguins parameters p^. and to zero. Comparison with data then 
gives an insight about the importance of nonperturbative charming penguin contributions 
and/or on the size of subleading terms as detailed below Eqs. ( 150]) and ( |5T1) . To reduce 
the hadronic uncertainties, the predictions for B — > XM branching ratios are normal- 
ized to dT(B — > Xrf/dEry. The predicted ratio of partial decay widths, Eq. (lT5"j) . de- 
pends on the light meson energy Em- In the endpoint region, however, the dependence on 



x M = 2E M /m B = 1 + 



m 



M 



Px)/ m l 



is a subleading effect. 2 We neglect this dependence 



and set xm = 1 in Tables [TTT] and HVl 

For the coefficients in the Gegenbauer polynomial expansion of the LCDAs 



(J)m(u, /i) = 6mm 



£<( W )C< 3 / 2 (2m 



n=l 



(49) 



we take the same values as Ref. except for the first coefficient in the Gegenbauer expan- 
sion of 4>k{x) for which we use the recent lattice QCD determination af (« = 2.0 GeV) = 
0.055 ±0.005 Q. Explicitly the remaining coefficients are (at « = 2 GeV): af = 0.23 ±0.23 
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K* 



a 



0.±0.15 36 



0.08 ± 0.13, [36|, al = 0.09 ± 0.15 [37 



and for lack of better information 



0.07 ± 0.08, af; 

il a = a\ and a% 



= 0.14 ±0.15, 
0.±0.2, with 



the higher coefficients in the expansion set to zero. 



2 For instance, the same p x cut corresponds to higher Em cut for heavier mesons. For p x < (2 GeV) 2 one 
has E n > 2.26 GeV for B -> ttX, while > 2.36 GeV for B — > 4>X (to be compared with m B o/2 = 2.64 
GeV). Thus m B /2 - E M ~ A with 1 - x M ~ 0(A/m B )- 
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MX 


Br(£?~ -> MX)/Bt(B -> 




A C p 




0.67 ±0.37 ±0.14 




-0.04 ±0.02 ±0.01 


CLU 


(4.1 ±2.1 ±2.6) x 10" 


-3 


0.64 ±0.10 ±0.10 




(1.0 ±0.5 ±0.3) x 10" 


-2 


-0.15 ±0.11 ±0.01 


P 


1.76 ±0.97 ±0.38 




-0.04 ±0.02 ±0.01 




(1.3 ±0.6 ±0.7) x 10" 


-2 


0.63 ±0.10 ±0.10 




(1.4 ±0.8 ±0.5) x 10" 


-2 


-0.17 ±0.11 ±0.03 


<t>*M 


(2.0 ± 1.1 ±0.1) x 10" 


-4 




»** 


(3.8 ± 1.8 ± 1.1) x 10" 


-3 


-0.72 ±0.13 ±0.20 



TABLE IV: Predictions for decay rates and direct CP asymmetries for charged B~ — » MX AS = 
semi- inclusive hadronic decays, which are the same as for corresponding B — > MX, B s — > MX 
given in Table HI1 The first errors are an estimate of the l/m& corrections, while the second errors 
are due to errors on the Gegenbauer coefficients in the expansion of the LCDA. 

Direct CP asymmetries, Eq. ( Jl~6l) . are nonzero only in the presence of nonzero strong 
phases. These can be generated nonperturbatively or by integrating out on-shell light quarks 
in a loop when matching full QCD to SCETj at NLO in a s . As in Ref. j^] we therefore use 
the NLO matching expressions for the Wilson coefficients Cf at fi = mi,, to have the leading 
contribution to the CP asymmetries, while performing the evolution to the hard-collinear 
scale /io ~ ^/Arrib at NLL. Note that this running cancels to a large extent in the ratios 
of the decay rates (only the running of a^(/i),n > 1 remains), giving in effect the Wilson 
coefficients with NLO accuracy at the hard-collinear scale /io [5|- We choose = 2 GeV for 
the perturbative predictions in Tables Hill and HVl 

The two errors quoted in Tables II III and IIVI are an estimate of subleading corrections 
and due to the errors on the Gegenbauer polynomial coefficients in the LCDA expansion 
( )49l) . Since the predictions are made to NLO in a s (rrib) but only to LO in 1 /nib, the largest 
corrections are expected to arise from the 1/m^ terms. These are estimated by independently 
varying the magnitudes of the leading terms proportional to \^ ct by 20% ~ 0(A/m b ) and 
the strong phase by 5°. This latter variation estimates the error on the strong phase arising 
from the uncalculated a s {m h ) /m b or ot 2 s {m h ) terms. A 100% error is assigned to predictions 
for branching ratios in color-suppressed tree and QCD penguin-dominated AS" = decays 
where the 1 /m& corrections are sizable compared to the leading results due to the hierarchy 
of Wilson coefficients. No prediction on CP asymmetries is given for these modes or for the 
QCD penguin-dominated AS = 1 decays. 

Next we confront the perturbative predictions with experimental data. Normalizing the 
BaBar results on semi-inclusive B — > KX branching ratios, Eq. ([1]), to Br (B — > X^) = 
(172 ± 21) x 10 -6 with the same photon energy cut E 1 > 2.34 GeV that was used for the 
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kaon momentum 381, one has 



T(B-/B° -> K-X) 1iq , nqn ,- m 

A_Z__J = o, 89 ±0 , 42 . (51 ) 

The central values of the measurements are substantially higher than the perturbative pre- 
dictions for B~ — > X _ X° S and _B~ — > if X^ modes, 3 given in Table IIHI They are still 
consistent within errors, with the discrepancies respectively at 3 a and 1.6 a levels for K + X 
and K°X modes, but do indicate that there might be substantial nonperturbative charm- 
ing penguin contributions (or very large l/m& corrections). Using isospin symmetry the 
charming penguin parameters, pf c and V^., in the two modes are the same. The three 
real parameters, \p^ c \, arg(p^), and V^,, can then be determined from the four observables: 
the two branching ratios (1501) . AST]) and the corresponding CP asymmetries once these are 
measured. This would also leave one observable as a consistency check. 

Since the CP asymmetries are not measured yet, this procedure is not possible at present 
without further approximations. Quite generally one expects that roughly \p^ c \ 2 ~ Vf c . As 
a starting point we thus take this relation to be exact and then extract \p£,\ as a function of 
arg(p^) from the K + X (K°X) decay width. This gives for the nonperturbative charming 
penguin to be about a factor 4 ± 3 (4.5 ± 2) larger then the perturbative prediction, i.e. the 
part of Q containing the C\ Wilson coefficient. 

Here the error is a sum of the experimental error and the error due to the variation of 
arg(p^) e [0,27r). Another way of presenting this result is through the ratio of nonpertur- 
bative and perturbative contributions to the decay width 



A c Pec 



h (s) 

a K 



2.2 ± 1.1 : K+X 
2.0 ± 1.5 : K°X, 



(52) 



that should be zero if the charming penguins are purely perturbative. The error on the 
ratios mostly reflects the variation due to a scan over the phase of p^ c . These values are 
very sensitive on the assumed relation between \pf c \ 2 and and should be taken as a 
rough guide only. Nevertheless, they show that there is experimentally a possible indication 
for sizable nonperturbative charming penguin. Two ingredients would help to clarify the 
situation significantly. First, the inclusion of chirally enhanced \ I Trip terms would show 



whether part of the discrepancy can be attributed to those terms [22|, [23j. Second, the cut 
on p*(K) should be lowered experimentally below the rather high value of 2.34 GeV used 
at present [lj|, so that one would be sure that the measurement is in the endpoint region, 
where our calculations are applicable, and away from the resonance region. 



3 The measurements are an average over charge and neutral B decays to K~ X (or K X) final state, but 
these are the same to the order we are working, see Table |U Decays to K X final state include also an 
incoherent sum with AS = decays into K a X~ a , which are, however, CKM suppressed and thus small, 
see Table El 
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Further tests are possible once more modes are measured. For instance, by using SU(3) 
flavor symmetry the charming penguin parameters for different modes can be related to only 
three real parameters, \p cc \, arg(p cc ) and V cc , making the framework even more predictive 
(at the expense of some accuracy due to SU(3) breaking). Note that in order to realize 
experimentally whether there are nonperturbative charming penguins and/or large l/m& 
corrections one does not need any symmetry arguments, just a comparison between our 
perturbative predictions and the experiment. To distinguish between the two sources of 
corrections, however, the flavor symmetries would most likely be needed. 

An interesting set of modes that could be used to settle the question concerning large 
1 / mi corrections versus nonperturbative charming penguin are the decays where no charming 
penguins are present. These are the AS = decay B — > (fiX and the color suppressed 
AS = 1 decays B — > uX, tt°X. They are experimentally more challenging since one would 
also need to measure the strangeness content of the inclusive jet. Namely, the related decays, 
the AS* = 1 decay B — > <pX and the AS = decays B — > ujX, tt°X, do contain charming 
penguins and should thus be distinguished experimentally from the first set of modes that 
does not contain charming penguins. Furthermore, the decays that do not receive charming 
penguin contributions might have substantial l/m^ corrections [5], so one may need to go 
one higher order in the expansion to have a definite understanding of experimental 

results. 

We next move to the decays involving rj and rf mesons. The expressions for the decay 
widths in the endpoint region have been derived in Section HVl with the final result given in 
Eq. fj48|) . At LO in l/m^, there appear two new shape functions that are specific to B — > 
rj^X decays. At lowest order in a S) S g (EM,^o) and S gg (E]\.j, /^o) defined in Eqs. ( )43l) . ( )46|) 
(or T g l and T gg in Eqs. ( 1B17|) and (1B20j) to all orders in a s ) describe gluonic contributions 
coming from one and two 0\ g operator insertions respectively, cf. Fig. [5l Little is known 
about these new nonperturbative functions, because of the lack of experimental data. At 
oresent only the B — > rj'X partial decay width with E v i > 2.218 GeV cut has been measured 
2, 0, S] • Normalizing to the B — > A7 decay width with the same E 1 cut gives 4 



AS = 1 and AS = 05^ v'X, t]X decay widths and direct CP asymmetries. Working at 
LO in a^VAmj) and neglecting the E M dependence of shape functions, one introduces only 
two new nonperturbative parameters specific to these decays, S g and S gg . With more data 
they could in principle be determined from experiment in the future. Note in particular that 
the charming penguin parameters entering the predictions can be fixed using SU(3) flavor 
symmetry from semi-inclusive decays into nonisosinglets. 

To have some future guidance on the size of the missing components, we list in Table 
IVl the purely perturbative predictions where we set both the charming penguin as well as 



4 If instead the same p v ' and p 7 cuts are used the decay widths ratio is 1.34 ±0.30. The difference compared 
to Eq. ([53")) reflects the effect of the large mass. 




In B 



77 



■OX decays there are 8 observables that are independent at LO in 1/m^: the 
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the new gluonic shape function contributions to zero and neglect NLO corrections from 
C>2g insertions. One striking aspect of that calculation is that the prediction for AS = 1 
B — > i]'X s decay falls remarkably short of the large measured value. Furthermore in this 
incomplete perturbative prediction the hierarchy between rjX s and rfX s decays is exactly 
opposite to the one found in two body B — > r]'K, r\K decays. While Br(£> — > rj'K) 3> 
Br(B — > rjK), the hierarchy in the incomplete prediction for the semi-inclusive decays is 
inverted. In the two body decays the hierarchy is well explained through the constructive 
and destructive interference of B — > rj q K and B — > rj s K contributions in B — > rj'K and 
B — > rjK amplitudes respectively due to 77 — rf mixing, if A(B — > rj q K) ~ A(B — > T] S K). 



This is exactly what is found in the limit of charming penguin dominance 2l|, |39|. In 
the semi-inclusive decay on the other hand there are no charming penguin contributions to 
B — > i] q X s at LO in l/m&, cf. Eq. ( 130]) . so that with large charming penguins there is no 
large hierarchy between B — > rj'X s and B — > r^X s decays. 

To be more quantitative, it is instructive to take three formal limits: (i) dominating 
charming penguins, (ii) dominating Oi g contributions, and (hi) the incomplete perturbative 
prediction with C 2 # — ► 0. If the amplitudes are dominated by charming penguins then (in 
the SU(3) symmetric limit) 

Br(fl- - ^ 2 



Bt(B- -> r/'X; 



tan 2 = 0.67. (54) 



Thus, as already argued above, in this limit there is no hierarchy between the two decays 
since Lipkin's argument of destructive and constructive interferences does not work for semi- 
inclusive decays. On the other hand, if Oi g contributions dominate then working at LO in 
a s (y/Arrib) we have 

Br(B~ -> rjX^) _ |(cos0/ % % - sin^JJ ® ( J + J) | 2 



Br(B- -> VX sS ) |(cos^^ + sin0/„>„) ® (J + |) 



|2 ' 



(55) 



Numerically this gives 1.2 x 10~ 4 if asymptotic LCDA are used, and 1.52 x 10~ 2 if SU(3) 
breaking is estimated by setting a^ s = instead. In the limit of dominant Oi g contribu- 
tions we thus have a similar large hierarchy between r]X s and i]'X s decays as in the two-body 
decays due to the destructive interference as apparent from Eq. (1551) . Finally, if the incom- 
plete perturbative calculation were a valid approximation, then we would have an inverted 
hierarchy between 7]X S and rfX s decays. This is due to a cancellation that is found between 
different terms in the B — > r(X s perturbative prediction. Due to small B — > r(X s decay 
width, however, this limit is phenomenologically excluded. 

In order to understand the relative size of charming penguin and the O^g contributions 
it is important to have a measurement of B — > i]X s decays. The relative size compared to 
the B — > i]'X s decay width is clearly different in the two extreme cases when only one of 
the two contributions is important. Comparing further with the other decays one should be 
able to determine all the nonperturbative parameters. As an exercise we set V Vs = |Pccl 2 
and take to be equal to obtained from B — > K~X, leading to a prediction for the 
normalized decay width Br(Z?~ — > rj' 'X~ u ) /Br(B — > X s ^) = 0.43 ± 0.25 with the variation 
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MX 


Bt(B~ -> MX)/Br(B -> 


Xsl) 


A C p 




(5.6 ±2.9 ±0.6) x 10" 


-2 


(-7.1 ±1.8 ±2.8) x 10~ 2 




(1.0 ±2.0 ±0.3) x 10" 


-2 


0.19 ±0.19 ±0.08 




(6.2 ±3.2 ±1.3) x 10" 


-2 


-0.38 ±0.10 ±0.10 


v'x da 


(2.4 ±1.2 ±0.7) x 10" 


-2 


-0.46 ±0.12 ±0.10 



TABLE V: Predictions for decay rates and direct CP asymmetries for charged B — > ryA , rj'X 
decays AS = 1(0) semi-inclusive hadronic decays given above (below) the horizontal line. The 
predictions equal also the corresponding B° — > r)X,r)'X, B° s — > r]X,r/X decay as given in Table HH 
The first errors are an estimate of the l/m^ corrections, while the second errors are due to errors 
on the Gegenbauer coefficients in the expansion of the LCDA. 

mainly due to the unknown strong phase of p^ c (while for the AS = decay we find 
Bt(B~ -» rfXj a )/Bi(B -> X,7) = 0.02 ± 0.02). This is still smaller then the measured 
value (}53l) . not surprising given the approximations made to arrive at it. Whether the 
difference is partially explained also by C>2g contributions should be clarified once more data 
are available. 

VI. CONCLUSIONS 

In the framework of SCET we considered semi-inclusive, hadronic decays B — > XM in 
the endpoint region, where the light meson M and the inclusive jet X with p 2 x ~ Amj are 
emitted back-to-back. This is an extension of the analysis done in Ref. \^ where we limited 
consideration to decays in which the spectator quark does not enter into the meson M. The 
contributions in which the spectator quark enters the meson M are power suppressed by 
1/ml in SCET. In this work we thus extend the SCET predictions at LO in 1/nif, to all decays 
where the spectator can enter either the jet X or the meson M. In SCET the four-quark 
operators factorize, which allows for a systematic theoretical treatment. After matching the 
full QCD effective weak Hamiltonian onto SCETi, the weak interaction four-quark operators 
factor into the heavy-to-light current and the n-collinear current. The forward scattering 
amplitude of the heavy-to-light currents leads to a convolution S of the jet function with the 
5-meson shape function, while the matrix element of n-collinear currents gives the LCDA 
for the meson M. The product of the two then gives the factorized form for the decay 
rates. The two nonperturbative functions, the convolution S and the LCDA, are the only 
nonperturbative inputs in the predictions for B — > XM decay rates at leading order in l/m*,. 
Furthermore, the same convolution S appears in B — > X s 7 decay and drops out in the ratio 
of B — > XM to the B — > X s 7 rate and in the predictions for direct CP asymmetries. Further 
work on higher order corrections would be useful in reducing the theoretical uncertainty. 

Nonperturbative charming penguin contributions can be included by the addition of one 
real and one complex parameter in the SU(3) symmetry limit. These parameters, which are 
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zero if the charming penguins are purely perturbative, can in principle be extracted from 
data. Thus by investigating decays without charming penguins, we can test whether the 
formalism is working. Then by looking at modes where the charming penguin can contribute, 
we can potentially see whether or not the charming penguin gives a large contribution to 
the decays. 

Decays where the light meson is an isosinglet rj or rf are special in that they receive 
additional contributions from gluonic operators. We consider these decays in detail, and 
show that the decay rate still factorizes, but there are two new shape functions which enter 
into the predicted rate. 

Using the available data, we performed an analysis of the semi-inclusive hadronic decays. 
While to date the data is limited, our preliminary analysis seems to indicate either large 
higher order corrections or a large contribution from non-perturbative charming penguins. 
With more data, it may be possible to distinguish the two possibilities as well as to extract 
the size of the nonperturbative charming penguin contributions. 
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APPENDIX A: TREE-LEVEL MATCHING 



The matching of the effective weak Hamiltonian in full QCD 

10 



G p 

71 



J2 A? (CiOf + C 2 Ol) - Al 9) (J2 CiOi + c g o g + c 7 o, 

j)=u,c i=3 



(Al) 



onto an SCETi one was calculated at NLO in a s (mb) first in Refs. 23J, and then in Ref. 19 



giving the SCETi effective weak Hamiltonian in Eq. (17!) . Here we list for reader's convenience 
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the tree level result of the matching 
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(A2) 



while the NLO results in our notation can be found in Appendix A of Ref. |5 



APPENDIX B: DERIVATION OF GLUONIC CONTRIBUTIONS 

In this appendix we provide details on the derivation of G(0) in Eq. (1351) and extend the 
results for (dT/dEM) g and (dT/dE^gg in Eqs. ( l4"5l) and (j47j) to all orders in a s (^f^Vm^ ) ). We 
start with the T-products 

g mcg) = (MX\i J d*xT{0 lg {Q),d§ {cg) {x)}\B), (Bl) 

that were already defined in Eq. (1521 and are also shown in Fig. El The subleading SCET 
Lagrangians appearing in Eq. (IB 1|) are 



=C (l# A^) ^ W< Pi ^n)^ + & PnWn) {Y}i ft?*) 



where the sum over light quark flavors q' is understood, and [4 



C S = Tr ( [^o, iD?] [iD^, WniDiM] }, (B3) 

with z£)^ = i£>^ + and W» = ^-P + P^ + ^-in ■ D us . 
For the calculation of it is useful to rewrite 

= £« a + =^g^-l^ig + ^ #^_±_!^ fl + • • • , (B4) 

with defined in Eq. ( 1541) . while the ellipses denote additional terms containing P± that 
do not contribute in our case. is then also split accordingly into 

QKMb) = (MX\i J d 4 xT{O lg (0)i£ { £a( b) (x)}\B)- (B5) 

In the SCETi to SCETn matching (where p 2 ~ Amj intermediate degrees of freedom are 
integrated out) we focus on the n fields in Q^ >a - The matching leads to two jet functions 
once hard-collinear modes are integrated out, and we obtain 

T{[B^] cd (0), [{qn) a {llPiiqn) b }u{x)} = i6(x_)8\x ± )± 
f dk 



mb (B6) 



-_ e - lk . X+ /2 [(TA^a^Aydj^ k ) + gab god fc _ } ] ^ + 

Z7T 
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with A;_ = n ■ k, x + = n ■ x. The ellipses are terms which do not contribute to rf' states. 
Tree-level matching gives for the jet functions Ji(u, kJ) = l/(Nk-) and J[(u, kJ) = 0. The 
J[ term does not contribute to G^, a since it leads to Tr(gA^) = 0. The remaining piece 
can be rearranged using color identities into 

g *f dk-dx + e _ lk _ x+/2 I du %L) 

47r J Air J u ( B7 ) 

x (MX\[q n ^P R Y^Y n gAi a (x)Yib v } [q nl ^ ± iq n ]u\B). 

The two terms in the square brackets are factorized in the sense that there are no soft 
gluon exchanges between the two terms - all the soft fields are in the first bracket. The 
communication between the two is only through the k_ and u convolutions with the jet 
function Ji(u, k_). 

Using the definition of the LCD A 



1 ft ab /til \ 

(M\(qn)mA\0) = -^EjufMMv) — ^), (B8) 



Sab 

to evaluate the matrix element from the second square bracket in (1B7[) we then have 



fe, = / duf M <p M (u) [ d ^ ± ^ ik - X+/2 -Hn 1 kJ) 

4tt J J Att u ( B9 ) 

(X\q n Y^Y n (0)^i(x + )P R YX(0)\B). 



x 



In simplifying the Dirac structure the identity 

eT ¥$]l^P R = ^IPr (BIO) 

was used, e 0123 = +1, and = e a ^ XcT n\n a /2. For the derivation of G^,b we notice that C$ b 
is a hermitian conjugate of C^ a - Using the hermitian conjugate of (IB6[) we finally have 



(Bll) 



x (X\q n Y^Y n (0)^gM(x + )P R Yib v (0)\B). 



Moving now to the calculation of Q cg , we first rewrite C c l g in a more useful form 

£« = ^Tr{[igB^JgB^][igB^,gAi v }} + --- , (B12) 

where the ellipses denote terms that do not contribute to Q cg . The matching from SCETi 
to SCET n gives 

T{i 9 m ± r(0),£$(x)} = K U 9(Ai u r d (x)5(x^)5\x ± ) / __e-^^/ 2 



-eg 
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X 

m b 



2tt 

(B13) 



duJ g (u, k-)e£" Tr [igB^igB^] u + 
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where again the ellipses denote terms that do not contribute for 7/ ) final states either because 
the collinear gluons are not in color singlet combination or they have incorrect parity. Using 
the definition of gluonic LCDA Eq. fl2TT) and the identity fIBlOj) . we then have 



'eg 



-WCi 



- u s m b 
4tc 



dufh& P (u) 



elk— 

47T 



-«-*+/2j fl ( U ,Jfe_) 



x (X\q n Y^Y n (0)^gM(x + )P R Y^b v (Q)\B). 



(B14) 



At tree level we have J g (u, k_ 

¥ P {u) -- 



- Jg(k_), independent of u. Since §p(u) is antisymmetric, 
-<&p(u), the matrix element Q cg vanishes at this order. The sum of the two 
contributions, in Eq. (IBllf) and Q cg in Eq. (IB140 then gives the result for Q(0) quoted 
in Eq. (155]). 

We next extend the results for (dT/dEM) g and (dT/dEM) gg in Eqs. (]4"5]) and f]4"T]) to all 
orders in q s (\/Sj). To do so, we redefine the heavy-to- light current J in Eq. ( |38l) to 
contain also the integration over hard momenta fractions u [since in general one may not be 
able to factor this dependence from the dependence on soft k_ momenta in the jet functions 

Ji, g (u,k-)], 



J M (0) 



du 



dk_dx + 

47T 



e- ik - x+/2 F M (k_,u)J H (0,x, 



(B15) 



The heavy current Jh is given in Eq. ( |33l) . while the hard-collinear kernel multiplied by 
LCDA, F M (k„,u), is given in Eq. ([36]). Unlike the current J in Eq. ([38]), the current J M 
in (1B15I) depends on the final state meson M through LCDA that are part of the F M (fc_, it) 
function. 

The derivation of the B — > XM decay width is now very similar to the one given in Section 
HVl Starting from the T product of heavy currents corresponding to one (D\ g insertion, 



Tg 1 {Em) = — / d?z(B\Tf H {z)J M (Si)\B), 
y rrib 



(B16) 



with Jh{z) 



' z (q n fiPiX^b^)(z) and J7" M (0) given in flB15j) . we use the factorization 



of n collinear quark fields from the rest at LO in l/m& to write 
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J P {1 + + m b - 2E M 
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x f g (l + ,r_)F M (r.,u) = 2^ I (E M , f i 



(B17) 

while the shape 



where the n-collinear jet function Jp(n + + it) was defined in Eq. (T4T 
function f g (l + ,rJ) was defined in Eq. ([42]) . Using the optical theorem we now have for the 
decay width contribution from single 0\ g insertion 



dT 



dE 



M / 9 



q2 



(B18) 



This extends Eq. ( 143]) to all orders in a s (^Am b ). We reiterate that the "shape" function 
J r g I (EM, yUo) n ow contains an integral over hard momenta fractions in the LCDA so that it 
depends on the meson M. 
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Defining similarly for the double 0\ g insertion 

T%{E M ) = — t d 4 z(B\TJ M \z)J M (0)\B), (B19) 



we have 



Disc.7^) = 2/ dMr . rfs .X m[ -> + + mi -2 BM + i e)] (B20) 
x ! duF M (r.,u) I dvF M ( s _,v)'f„(l + ,r_, s _) = 2m i r"(E u ,^), 



where the shape function f gg (l+, r_, s_) that depends on three soft momenta was defined in 
Eq. (I4"5l) . For the double 0\ g insertion contribution to the decay width we then have 

{^A =^rn\x M Tf g {E M ^,)\\fC lg \\ (B21) 

\dE M Jgg 27T yy 1 y| 

which extends Eq. ( I4T1) to all orders in a^Amj), with the "shape" function jF^(S M ,/i ) 
again depending on the meson M through the LCDA. 
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